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The temporal response of two-wave-mixing in photorefractive InP:Fe under a dc electric field at 
different temperatures has been studied. In particular, the temperature dependence of the charac- 
teristic time constant has been studied both theoretically and experimentally, showing a strongly 
decreasing time constant with increasing temperature. 
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I. INTRODUCTION 

The photorefractive effect leads to a variety of nonlin- 
ear optical phenomena in certain types of crystals. The 
basic mechanism of the effect is the excitation and re- 
distribution of charge carriers inside a crystal as a result 
of non-uniform illumination. The redistributed charges 
give rise to a non-uniform internal electric field and thus 
to spatial variations in the refractive index of the crystal 
through the Pockels effect. Significant nonlinearity can 
be induced by relatively weak (|iW) laser radiation. Phe- 
nomena such as self-focusing, energy coupling between 
two coherent laser beams, self-pumped phase conjuga- 
tion, chaos, pattern formation and spatial soliton have 
attracted much attention in the past 20 years 0, Q . 

Among photorefractive crystals, semiconductor mate- 
rials have attractive properties for applications in optical 
telecommunications such as optical switching and rout- 
ing. This is due to the fact that they are sensitive in 
the infrared region and their response time can be fast 

Two- wave-mixing is an excellent tool to characterize the 
photorefractive effect in these materials 0, la, Lj by de- 
termining the gain of amplification under the influence 
of the applied field, impurity densities, or grating period. 
Some semiconductors, like InP:Fe, exhibit an intensity 
dependant resonance at stabilized temperatures [HQ. 

In this paper, we analyze the temperature dependance 
of Two- Wave-Mixing (TWM) characteristic time con- 
stant, theoretically at first and eventually against exper- 
imental results. We propose a formal description of the 
temporal evolution of carrier densities in the medium, 
linking them to the TWM gain temporal evolution. 



II. TIME DEPENDANT SPACE-CHARGE 
FIELD IN InP : Fe 

The basic principles of the photorefractive effect in 
InP:Fe are well known |(|- It involves three steps: pho- 



toexcitation of trapped carriers into excited states, mi- 
gration of excited carriers preferentially towards non- 
illuminated regions and capture into empty deep centers. 
This leads to the formation of a local space-charge field 
E sc and thus to the modulation of the refractive index. 
The modulated refraction index is then able to interact 
with the beams that have created it. When the modula- 
tion stems from beam interference as in two wave mixing, 
an energy transfer between beams may occur. 

The principle of two-wave-mixing is to propagate si- 
multaneously in a photorefractive crystal two coherent 
beams, which have an angle between their directions of 
propagation. This phenomena is governed by the follow- 
ing system of coupled nonlinear differential equations: 
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where I p is the pump intensity, I s is the signal intensity, 
and Iq is the total intensity equal to the sum I s 



a is the absorption coefficient (assumed here to be the 
same for pump and signal). In a photorefractive crystal, 
r takes the following form [5|: 
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where n is the refractive index, r e ff is the effective 
electro-optic coefficient, A is the beam wavelength in vac- 
uum and Im{E sc } is the imaginary part of the space- 
charge field ii^c (the ? shifted component of E sc with re- 
spect to the illumination grating). The expression of E sc 
will derived in the following lines. 9 is the angle between 
the two beams. 

In order to evaluate the photorefractive gain Tq given 
by equation © , the space-charge field E sc has to be cal- 
culated from the modified Kukhtarev model pj, taking 
into account both electrons and holes as charge carriers. 
We chose a model with one deep center donor, two types 
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of carriers (electrons and holes) Q , considering variations 
only in one transversal dimension (x) as described by the 
following set of equations: 
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Figure 1: Temporal evolution of electron (a) and hole den- 
sities (b) under uniform illumination for A = 1.06 |J.m, Io = 
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where E is the electric field, n e and ph are the elec- 
tron and hole densities in the respective conduction and 
valence bands, tit = Fe 2+ is the density of ionized oc- 
cupied traps, pt = Fe 3+ is the density of neutral un- 
occupied traps, J n and J p are respectively the electron 
and hole currents. Nt, and Na are respectively the 
densities of iron atoms, the shallow donors and the shal- 
low acceptors. The charge mobilities are given by \i n 
for electrons and \x p for holes, the electron and hole re- 
combination rate are respectively c„ and c p , T is the 
temperature and kb is the Boltzmann constant. The di- 
electric permittivity is given by e while e is the charge 
of the electron. V app is the voltage applied externally to 
the crystal of width d. The electron and hole emission 
parameters are e„ and e p depend on both thermal and 
optical emission as described by: 
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where the thermal contribution to the emission rate co- 
efficient is e th and the optical cross section of the carriers 
is given by a, I(x) is the spatially dependent intensity of 
light due to the interferences between pump and signal 
beams and hv is the photon energy. 

For sufficiently small modulation depth m, intensity 
and all carriers densities may be expanded into Fourier 
series interrupted after the first term : 



A(x) = Ao + Axe 



(7) 



where A(x) takes the role of /, n e , ph, nT, Pt and K g 
the spatial frequency of the interference pattern. So the 
light intensity can be written for the average intensity Iq 
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In the following, we have calculated the temporal evo- 
lution of carriers density and we will look forward to find- 
ing the temporal evolution of the space charge field under 
these hypothesis, i.e. considering only the zero'th and 
first order of the Fourier expansion. 

The zero'th order corresponds to an uniform illumina- 
tion ( I(x) = Iq). The space charge field is thus equal to 
zero and the local field is uniform and equals the applied 
field E, 

state are known to be equal to 

tively 0, where ut , Pt are the density of occupied and 
unoccupied traps at steady state. 

The electrons and holes densities in transient regime 
when an uniform illumination is established, are calcu- 
lated by solving equations i|4d|) and (|4*e|l , assuming that 
Iq = for t < and at time t = 0, the carriers density 
values are equal to hq and po at thermal equilibrium, 
without any optical excitations ■ We obtained the fol- 
lowing solution: 
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The temporal evolution of carrier densities under uni- 
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form illumination is illustrated in figure ^ Our model 
confirms the fact that the carrier densities evolution 
grows exponentially. The rise time of carriers genera- 
tion is on the order of nanosecond time scale for a beam 
intensity of a few mW per cm 2 . Without presence of 
the beam, the electron density is greater than the hole 
density because electrons are mostly generated thermally 
while holes are generated optically |5j. 

For a modulated intensity (first Fourier order), by us- 
ing the set of equations @, the space-charge field can be 
approximatively expressed at steady state as [jj: 



and 



The subscript indexes n and p are related to the elec- 
tron and hole contributions respectively. r„ and w n are 
given by: 
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where E Q , E d = K g ^f and E q = ^-^^ are the 
space charge field under uniform illumination, the diffu- 
sion field and charge-limiting field respectively. I res — 

e th n 

— " ' T P is the resonance intensity defined as the inten- 
a p .p To .h.v •> 

sity at which holes and electrons are generated at the 
same rate. 

From equation JUJ, we observe that the space charge 
field is purely imaginary, when the illumination Iq equals 
I res - Above resonance, the hole density is higher than 
the electron one, mainly because the holes cross section 
is stronger than for electrons. The result is that charge 
transfer mainly occurs between iron level and the valence 
band. Below resonance, when electrons are dominant, 
the iron mainly interacts with the electrons and the con- 
duction band. 

In transient regime for a modulated intensity, the dy- 
namics of the space charge field is calculated by consid- 
ering an adiabatic approximations ,9J , a concentration of 
electrons and holes densities reaches instantaneously the 
equilibrium value which depends on the actual concentra- 
tion of filled and empty traps, so we set: ^jrp = = 0. 
We assume that the electrons are excited thermally while 
the holes optically Q . In the low modulation approxima- 
tion, some algebraic manipulations of the set of equations 
(gj lead to: 
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where Emu = — is the mobility field, Tdi, n is the 
electron dielectric relaxation time depending on intensity 
and temperature which can be written as: 
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is the thermal parameter equal to: 
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where — is the effective masse of electron, E na is the 
apparent activation energy of the electron trap, cr£° is the 
electron capture cross section. The value of this param- 
eters are determined experimentally jlCj . 
To obtain t p , w p and Tdi, P , the index p should be sub- 
stituted for the index n in the equations l|14f> . i|15|) and 

m- 

From equations (|12|l to 117|L it is possible to deduce, as 
was done previously 0, that the time constant r g is 
real if electron emission is equal to the hole emission. 
That is, in the case of InP:Fe, electron thermal-emission 
is equal to the holes optical-emission. This allows to in- 
fer a link between the behavior of InP:Fe as a function 
of intensity and temperature. 

It will be the aim of the following sections to confirm this 
link, both theoretically and experimentally 



E 1 (t)=mE sc [l-exp(- — )] (12) 

where r g is a complex time constant, which can be 
rewritten by separating its real and imaginary parts. 
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III. GAIN DYNAMICS 

The photorefractive gain T is the main parameter that 
can be determined by two- wave-mixing. It quantifies an 
energy transfer from the pump beam to the signal beam 
and is proportional to the imaginary part of the space- 
charge field. 

The gain value depends on different parameters like ap- 
plied electric field, iron density Nt, pump intensity and 
temperature |ll| . Our work concentrate on the study of 
the gain dynamics versus temperature and we particu- 
larly analyze the dependance of the rise time on temper- 
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Figure 2: Temporal evolution of local gain at different inten- 
sities at T = 297.fr: (e^ h = lG.Sls" 1 calculated from equa- 
tion UTTl . A = 5 \im, Eq = 10 kV/cm, crystal thick- 
ness L = 12mm : (a) Iq = 25.5 mW/cm 2 ~ Ires, (b) 
To — 15 mW/cm 2 , (c)7o = 50 mW/cm 2 . r r : characteris- 
tic time constant of amplification. 



ature. 

The stationary value of the photorefractive gain at dif- 
ferent temperatures is given by equation J3J) where E sc 
is given by equation (|llfl . This expression shows that a 
maximum gain is reached when Iq = I res . This maxi- 
mum corresponds to an intensity resonance p|- 

We studied theoretically the temporal gain behavior 
using the standard definition given by equation i|18[) de- 
duced from equation i|12|) by developing E sc and r g . 

r = r i + exp — x L U ] is 

r sin ip 

where 

ip: argument of stationary space-charge field (E sc = 
\E sc \expiip) 

t: amplification's characteristic rise time. 

Our theoretical simulations produce the curves repre- 
sented on figure illustrating the evolution of photore- 
fractive gain as function of time for three different pump 
beam intensities: at resonance, below and above reso- 
nance for the same parameters as in figure ^ We see 
that the gain amplitude differs from each intensity to an- 
other, it takes the maximum value around resonance. 

As a next step, we studied theoretically the TWM gain 
time response as a function of temperature. For an easier 
comparison with experimental results, in the following, 
the response time is r r will be considered as the time 
interval necessary for the gain to reach 90% of the first 
maximum of each curves as shown in figure [21 The re- 
sponse time r r versus temperature are given in figure |3| 
for three distinct intensities, along with a fourth fitted 
curve showing the time constant at resonance intensities. 
We observe that the response time quickly decays as tem- 
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Figure 3: Characteristic time of amplification versus temper- 
ature at three different intensities. The time constant at the 
resonance intensity is given for each temperature by the dot- 
ted line. 
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Figure 4: Local (a) and average gain (b) versus time, Jo = 
25.5 mW/cm 2 , ef = 16.31s" 1 , A = 5 urn, E = 10 kV/cm. 

perature increases. At resonance, r r is larger because the 
space charge field is high and it consequently necessitates 
more charge to accumulate. 



IV. AVERAGE GAIN 

The theoretical curve shown in figure |2 illustrates the 
temporal evolution of the local gain. For the InP:Fe sam- 
ple, the absorption coefficient at A = 1.06/im being ap- 
proximately equal to 1cm -1 . Owing to this absorption, 
the mean intensity decrease along the z axis propagation. 
The exponential gain would result from an integration 
over the optical thickness, as described in equation 1191 

Y = y ( T(z).dz (19) 
L Jo 

The figure 0] shows temporal evolution of local and av- 
erage gain for L = 12mm crystal thickness for the same 
intensity; the average gain is lower because the intensity 
absorption is taken into account. 
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• 10=28 niw/cm 2 

• 10—62 mw/cm 2 

• 10=97 mw/cni 2 




T(K) 



Figure 5: Average gain characteristic time versus temperature 
at three different intensities. The time constant at the reso- 
nance intensity is given for each temperature by the dotted 
line. 




Figure 6: Two-wave-mixing configuration. 



Because of the absorption the resonance intensity for 
average gain is higher than the local one for the same 
temperature. As for the local gain, we calculated numer- 
ically the characteristic time, in the same way. The re- 
sults are shown in figure^] We compare the results in fig- 
urelHl we observe the following differences: the resonance 
peaks are slightly widened because is reached within the 
example for various input intensities and the peaks are 
shifted towards high intensities again because of absorp- 
tion. These conclusions can arise from figures |3] and [S] 
although they show the rise time as a function of tem- 
perature. Indeed, our calculations show that the pho- 
torefractive gain and rise time are linked, so that the rise 
time is the slowest for the highest gain (i.e. at resonance); 
since more charges need to be accumulated. 



EXPERIMENTAL VALIDATION 




Figure 7: Gain dynamics at 288K at different intensities: 
(a)62 mW/cm 2 ,(b)97 mW/cm 2 (c)28 mW/cm 2 . 



We perform standard two- wave mixing experiments in 
co-directional configuration as shown in figure [5] Pump 
and signal beam intensities ratio is set to (3 = 50 and the 
angle between pump and signal is 29 = 12° corresponding 
for an space grating A = 5/im. The experiments are 
performed with a CW I.06(xm YAG laser. 

An electric field (lOkV/cm) is applied between the < 
001 > faces of InP:Fe crystal (5 x 5 x 12mm 3 ). The light 
beam is linearly polarized along the < 110 > direction 
and propagates along the < 110 > direction (12mm). 
The absorption constant as measured by spectrometer 
is close to 1cm -1 at 1.06//m. Crystal temperature is 
stabilized by a Peltier cooler. 

Transient behavior is analyzed by measuring r r as was 
done in figures and 03 Figure shows the results ob- 
tained for three different intensities from one side to the 
other of the resonance (the oscillations seen on figure [7| 
are attributed to the experimental noise and the curves 
are assumed to correspond to the first order responses). 
Experimental results concerning the TWM time constant 
are given on figure[S] For high temperature, r r decreases 



for all intensities. 

Note that, for both theory and experiments, t t value 
decreases from 300 to 50 ms for an increase temperature 
of 10°C-showing a good quantitative result. The dis- 
crepancy observed between figure|S]and[3]is partially cor- 
rected by taking into account the gain integration along 
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Figure 8: The r r versus temperature at different intensities. 



6 



the beam path inside the crystal, as shown in figure HJ 
showing a widening of the curves. We attribute the differ- 
ence observed in terms of gain maximum value to lack of 
precision in the knowledge of the crystal's physical con- 
stants such as photo-excitation, cross section and dopant 
concentration. 

VI. CONCLUSION 

We have studied the dynamics of TWM in InP:Fe as 
a function of intensity and temperature. A theoretical 



analysis shows that the gain coefficient oscillates when 
an intensity lower or higher than the resonance intensity 
is used. At resonance the gain grows exponentially. 

The experimental study shows that the crystal absorp- 
tion prevents the oscillating behavior. We have shown 
that the gain rise time is strongly temperature depen- 
dent. Experimentally the gain rise time is 10 times 
shorter at 295K than at 285K for low intensities. 

According to experimental and applications needs, the 
temperature as well as the intensity can be used to tune 
the photorefractive response time. 
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